Spacing Characterization in Al-Cu Alloys Directionally Solidified Under Transient 

Growth Conditions 
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We study spacing selection in directional solidification of Al-Cu alloys under transient growth 
conditions. New experimental results are presented which reveal that dendritic spacing versus solid- 
ification rate evolves in an almost step-wise fashion, consistent with previous theoretical predictions 
of Langer and co-workers. Phase field simulations of directional solidification with dynamical growth 
conditions approximating those in the experiments confirm this behavior. Changes in dendrite arm 
spacing is shown to be consistent with dendrite instabilities confined, initially, to sub-domains, rather 
than the entire system. This is due to the rapid variation in growth conditions, which prevent the 
system from adapting as a whole but, rather, in a succession of quasi- isolated domains. 
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I. INTRODUCTION 

Solidification microstructure is the starting point of 
any casting operation. Dendritic spacing and morphol- 
ogy established during casting often sets the scales of 
the downstream microstructure during manufacturing of 
alloys. This is particularly true in emerging technolo- 
gies such as twin belt casting, where a reduced amount 
of thermomcchanical downstream processing reduces the 
possibility of modifying microstructure length scales from 
that determined at the time of solidification. 

Predicting columnar microstructure in cast alloys has 
been traditionally studied in the context of Bridgeman 
type directional solidification conditions. Most studies 
have focused on the problem of primary and secondary 
arm spacing in dendrite arrays of thin liquid films of or- 
ganic alloys, directionally solidified under steady-state 
cooling conditions, i.e. a sample is pulled through a 
constant thermal gradient at a constant pulling speed. 
Careful experiments on steady state directional solidifi- 
cation reveal a reproducible correlation between spacing 
and pulling speed [l|, 0). Studies of steady state direc- 
tional solidification have developed so-called geometric 
models to relate spacing to solidification processing pa- 
rameters such as the pulling speed V, thermal gradient 
G and the alloy concentration c Q . 

In geometric theories of spacing selection the structure 
and mathematical form of the dendrite arms is first as- 
sumed and a subsequent consistency relation is derived 
for the arm spacing (also referred to as "wavelength" in 
the literature). The construction and assumptions that 
went into setting up the geometry of the dendrite array 
lead to at least one phenomenological parameter that is 
then fit to match the theory onto specific experimental 
spacing selection data [3|-l5|- While useful in helping to 
elucidate some aspects of spacing selection, such theo- 
ries lack the fundamental element of microstructure pre- 
dictability: the ability to self-consistently generate the 
morphology of the structure they are trying to predict 
anything about. It is also not clear how such theories 



hold up to a change of conditions away from those of the 
experiments they were constructed to model. 

While steady state directional solidification is an im- 
portant academic paradigm, it is not a realistic repre- 
sentation of the conditions prevalent during industrial 
casting, which normally occurs under rapidly chang- 
ing growth conditions. Near a chill surface or, indeed, 
throughout the entire sample thickness in the case of very 
thin strips, the thermal gradient and solidification speed 
are neither constant nor independent of each other. Rec- 
ognizing this, several research groups have attempted to 
extend geometrical models to include unsteady state pro- 
cessing conditions. For example, Garcia and co-workers 
Q and Kirkaldy and co-workers [ij have attempted to 
link the spacing of dendrite arms to the cooling rate q. 
Once again, as with any geometric theory, phenomeno- 
logical parameters are introduced to fit the model with 
experiments, although in this case the fits are not as good 
as in the steady state case. 

The importance of transient thermal history process- 
ing conditions in establishing as-cast microstructures has 
been apparent since the theoretical work of Warren and 
Langer M, |8| and the experimental studies of Losert and 
Huang 0, on alloys of succinonitrile (SCN). War- 
ren and Langer performed an analysis on the stability of 
dendritic arrays 0, Q and showed that they remain sta- 
ble to period doubling instability over a range of pulling 
speeds, contrary to the predictions of any geometric the- 
ory, steady state or transient. Their work suggested a 
band of available spacings rather than a unique selec- 
tion. In turn, Huang et al @ showed that by chang- 
ing the rate of the pulling speed it is possible to obtain 
different spacings for a given set of final growth condi- 
tions. Moreover, Losert et al observed that under 
a more gradual change in pulling speed the spacing still 
remained stable over a range of pulling speeds, consistent 
with the Warren and Lager predictions. They also noted 
the presence of a sharp increase in spacing which, con- 
sistent with the predictions of Warren and Langer, could 
be attributed to period doubling. 
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The results in the literature seem to point to two ex- 
tremes. Under steady state conditions, i.e. dV/dt = 0, 
the dendrite arm spacing appears to scale as a power- 
law of the pulling speed V, a result at least borne out 
qualitatively by geometrical models. On the other hand, 
under transient solidification conditions, dendrite spac- 
ing and structure seem to depend strongly on transient 
history and initial conditions, at least in the idealized 
setting of a linear stability analysis or for well controlled 
experimental SCN dendrite arrays. The lack of unified 
theory to explain both these regimes likely points to an 
incomplete picture of a fundamental physics underlying 
microstructure selection in solidification. It also points 
to the need for a robust theory and modeling formalism 
that predicts the evolution of dendritic morphologies and 
growth rate self-consistently, as function of only the in- 
put material parameters and cooling conditions, steady 
state or transient. 

Phase field theory has emerged in recent years as 
promising candidate of a fundamental and self-consistent 
theory for modeling solidification microstructures. The 
first simulations to test spacing versus pulling speed in al- 
loys date back to the work of Warren and Boettinger [ll[ , 
who found a monotonic band of spacings versus pulling 
speed. The small system size used, however, precluded 
a quantitative comparison with experiments. Nowadays 
the phase field methodology has become more quantita- 
tive by "marrying" simulations of phase field models in 
the so-called thin interface limit [l2l . IT3I ] with novel sim- 
ulation techniques like adaptive mesh refinement [21] . A 
first step using phase field models to quantitatively model 
spacing in directional solidification was taken by Green- 
wood et. al [H[ in 2D and Dantzig and co-workers in 3D 
[l6j . These works modeled steady state directional solid- 
ification in SCN alloys and found very good agreement 
(in the 2D limit) with the 2D steady state spacing ex- 
periments. These studies suggested that, at least under 
steady state (i.e. Bridgeman growth) conditions and one 
type of initial condition (morphologically noisy initial in- 
terface) , there could be a single crossover scaling function 
interpolating between the two power-law spacing regimes 
seen experimentally and modeled semi-empirically by ge- 
ometrical models. 

Despite the success of phase field modeling in predict- 
ing steady state spacing, as well as other steady state 
properties such as cell tip structure [17], the methodol- 
ogy has not been used systematically to explore spac- 
ing under transient solidification conditions. Indeed the 
ability to model cell, dendrite and seaweed structure, ki- 
netic and surface tension anisotropy, different mobility, 
different thermal conditions and different initial condi- 
tion makes phase field modeling an ideal theoretical test 
ground to explore transient spacing development and how 
it may relate to the steady structures. 

This paper reports new experiments that study pri- 
mary spacing selection in directionally solidified Al-Cu 
alloys cooled under transient conditions closely related 
to those encountered in strip casting of Al alloys. The 



transient thermal gradient and interface speed are mea- 
sured and correlated to measured dendrite spacing. Our 
results are shown to be inconsistent with steady state or 
transient geometric theories. Instead they are more con- 
sistent with stable ranges of spacings versus front speed, 
connected by rapid changes in spacing at particular in- 
terface velocities. We also present new two dimensional 
phase field simulations that exhibit the same behavior as 
our experimental data. Analysis of the phase field sim- 
ulations is used to shed some light on the morphological 
development of dendrite arms during solidification under 
transient conditions. In order to manage the length of 
this paper, a sequel paper to this one will further explore 
the theoretical implications of the phase fields simula- 
tions presented here, focusing in particular, on the the- 
oretical connection between transient behavior reported 
in this paper to the steady state behavior previously re- 
ported. 



II. EXPERIMENTAL PROCEDURE 

As-received Al-0.34wt%Cu samples were used to study 
solidification microstructure evolution under transient 
cooling conditions. The experimental set-up is shown 
in figure ([T]). It consists of a cylindrical stainless steel 
crucible that is water jet cooled from below to promote 
upward solidification. The crucible is shielded by a cylin- 
drical alumina insulation with a thickness of about 10mm 
to prevent radial heat extraction. The inner, outer and 
bottom parts of the crucible were covered with a thin 
layer of sprayed boron nitride in order to reduce the heat 
extraction through the walls as well as providing a more 
uniform chilling surface at the bottom. The pressure of 
the spray was chosen high to prevents the formation of 
bubbles at the water-chilling wall interface due to the 
local vaporization. 

The temperature is measured at different heights from 
the bottom with K-type nickel-chromium based sheathed 
thermocouples that are tightened along a plate and are 
inserted into the melt through the top of the crucible as 
shown in the figure. The diameter of the chilling surface 
is set to 50mm (i.e. the bottom plate) and the thickness 
of the chilling surface to 3.5mm. The large chilling diam- 
eter helps to reduce the influence of the thermocouples' 
diameter, about 1mm, on the solidification process. A 
set of thermocouples were aligned in the axial (vertical) 
direction starting at 1mm away from the chilling surface 
and separated from each other by 1mm. In order to as- 
certain the one dimensionality of the heat flow in the 
vertical direction, an additional thermocouple was posi- 
tioned 12mm radially from the aligned thermocouples. 

The output from the thermocouples is acquired 
through a NI SCXI-1600 data logger and the calibration 
of which was set at the melting point of pure aluminum 
and pure zinc. Before pouring the melt into the crucible, 
the crucible along with the surrounding insulation and 
the alloy were heated in the same furnace to a temper- 
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ature of 1.1 times the liquidus temperature of the alloy, 
to compensate for the heat loss during the experiment. 
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FIG. 1: Schematic of the set-up used to obtain upward direc- 
tional solidification. 



The cooling curves obtained from the thermocouples 
define a region of radially uniform temperature that 
varies essentially in the growth direction. This region 
is referred to as our region of interest. Outside this re- 
gion, there is a temperature gradient towards the wall 
such that the unidirectional heat flow breaks down. Far 
from the bottom chilling plate, outside the region of in- 
terest, the solidification microstructure mostly consists 
of equixed dendrites. Within our region of interest the 
solidification microstructure is columnar in nature (i.e. 
oriented cells, dendrites or seaweed branches). We dis- 
regard any data outside our region of interest. As will 
be shown below, the grains examined within the region 
of interest are large enough in the direction transverse to 
the heat flow to disregard boundary effects. Only data 
from grains with transverse size larger than 1mm are re- 
ported here. In each sample solidified, three different 
directionally solidified grains emanating from the chill 
surface were analyzed for reproducibility and to provide 
statistical error bars we report in our results. 

Figure ([2]) shows the longitudinal morphology of the 
dendrite microstructure cut out from a grain in one 
of our samples. The microstructure appears cellular 
in nature, while there is some evidence of seaweed-like 
tip-splitting throughout the image. The microstructural 
length scale in the direction transverse to that of the 
heat flow (right to left in the figure) was analyzed at 
different distances from the chilling surface using power 



spectral analysis. A typical power spectrum is shown in 
the bottom frame of figure ©. The x-axis of the power 
spectrum denotes the frequency of the corresponding 
wave vector. The main peak is associated to the primary 
spacing and is consistent with would be obtained by 
the ASTM line intersection method. The smaller, high 
frequency, peaks account for smaller wavelengths such 
as sidebranches, splitting tips or the discreteness of 
the pixels in an image. The smaller frequency peaks 
correspond potentially to longer-wavelength interactions 
between dendrites branches. Figure ([3]) shows a typical 
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FIG. 2: (Top) cutaway of the longitudinal solidification mi- 
crostructure, growing from left to right. (Bottom) Unidimen- 
sional power spectrum of a transverse cut extracted from the 
longitudinal microstructure at a position 13mm from the chill- 
ing surface. 



image of the corresponding transverse microstructure, 
cut away from a grain at distance 15mm from the 
chilling surface. The bottom frame of figure <j3j> shows 
a polar plot showing the average length scale versus 
angle in the data in the top frame, obtained from a 
digitization of the corresponding cross section in the top 
frame. The polar plot was obtained by averaging the 
two-point direct correlation function from data such as 
that in the top frame of figure ([3]) along lines of different 
orientation. The polar -or "rose"- plot in figure @ 
represent the mean inter-branch spacing in the data 
along different orientations. Details of the procedure by 
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FIG. 3: (Top) Typical cross section of transverse dendrite mi- 
crostructure. (Bottom) polar plot representing average length 
scale versus orientation in the transverse image in top frame. 



which the polar plot is constructed can be found in [1 



III. PHASE FIELD MODELING 

Phase field simulations modeled the Al-Cu alloy system 
in the dilute limit of the phase diagram, which comprises 
straight solidus and liquidus lines of slopes of m/k and 
m, respectively. The equilibrium interface concentration 
jump at any temperature is thus given by partition rela- 
tion c s = kci , where c s (c; ) is the molar concentration of 
impurities at the solid (liquid) side of the interface, and 
k is the partition coefficient. 

Simulations formally neglected the latent heat by im- 
posing the temperature field by the form T(z, t) = Tq + 



G{t){z -zo- f*V p {t')dt r ), where T(z o ,0) = T , a ref- 
erence temperature, while G(t) and V p (t) are the local 
thermal gradient and pulling speed, respectively. These 
were extracted from our experiments as discussed fur- 
ther below. We neglected convection in the liquid, as 
the transport of impurities becomes governed essentially 
entirely by diffusion Moreover, since diffusion of 

impurities in the solid is typically several orders of mag- 
nitude lower than in the liquid, diffusion in the solids is 
neglected. 

Under the above assumptions, mass conservation 
across the interface takes the form cj(l— k)v n = — Dd n c\i, 
where D is the solute diffusion coefficient in the liquid 
and d n \i is the partial derivative in the direction nor- 
mal to the interface, taken on the liquid side. The tem- 
perature at the interface, which is assumed to be in lo- 
cal equilibrium, is given by the Gibbs- Thomson relation 
T = T m — \m\ci — Tk — v n /fik, where T m is the melt- 
ing temperature of the pure material, T = jT m /L is the 
Gibbs-Thomson coefficient, 7 is the interfacial free en- 
ergy, L the latent heat of fusion per volume, k is the 
interface curvature, v n is the normal interface velocity 
and Hk is the atomic mobility at the interface. 

The underlying crystalline structure that defines the 
anisotropy of surface tension (or interface mobility) is 
modeled by through a commonly used fourfold symme- 
try anisotropy function a(h) = 1 — 3e + 4e(n^ + fly + nf), 
where e is the anisotropy strength and n is the unit nor- 
mal at the the interface. In two dimensions this function 
can be re- written as a(6) = 1 — ecos(4#), where 9 is the 
angle between the normal direction to the interface and 
an underlying crystalline axis (e.g < 100 > in a cubic 
crystal). 

Taking as reference the concentration on the liquid side 
of the interface c° = co/k (where Co is the average alloy 
concentration), the following standard one-sided sharp 
interface directional solidification kinetics are modeled: 

d tC = DV 2 c (1) 
d(l-k)v n = -Dd n \+ (2) 
c;/cf = 1 — (1 — k)ndoa(n) 

-Q.-k)(z-J V p (t')dA /l T - (1 - k)Pv n (3) 

where do = T/ATq is the solutal capillarity length, 
AT = \m\(l - k)c° the freezing range, l T = AT /G 
the thermal length, and j3 — l/(/ifcATrj) the kinetic coef- 
ficient. 

The phase- field model employed to emulate Eqs. ([TJ- 
§3§ is designed for quantitative simulations through the 
use of a thin interface analysis developed by Karma and 
co-workers [TU, [l3[ . This analysis makes it possible to em- 
ulate a specified capillary length and kinetic coefficient to 
second order accuracy in the ratio of the interface width 
to capillary length. The details of this model have been 
presented and discussed at length in Refs. [TH, [l3| and 
thus only a brief description is included below. A gen- 
eral review of the phase-field method can be found in 



5 



MM- 

A scalar phase field parameter <f> is employed, which 
takes on a constant value in each phase and varies sharply 
but smoothly across a diffuse interface. The phase field is 
used to interpolate the free energy density and mobility 
between the bulk phases. Its equation of motion guaran- 
tees the system evolves towards a minimum of the free 
energy of the system. We define a phase-field variable 
which takes the value (f> = l((f> = — 1) in the solid (liquid). 
The concentration c(x, t) is characterized by through a 
generalization of the field U — (c — c°)/(c°(l — k)), which 
represents the local supersaturation with respect to the 
point (c°,Xb), measured in units of the equilibrium con- 
centration gap at that temperature. This generalized su- 
persaturation field is given by 



U = 



1 



c/c° 



1 - fc V (1 - 0)/2 + fc(l + 0)/2 



- I 



(4) 



In term of the fields c, (f> and U, the phase- field model 
referred to above is given by 



+0-0 3 - A(l -(f> 2 ) 2 (U + 
fl + k 1-k ,\8U 



h 



(5) 



2 = V 



aw (l + (I - k)U)h 



dt 



q{4>)DVU (6) 

1 + {1 - k)U \d<f> 
2 ) ~dt 



where Zj nt = J Q V p dt' is the interface position, h = 

— (V0)/(| V0|) defines the unit vector normal to the in- 
terface, r(n) = To -a 2 (n) is the phase-field orientation de- 
pendent relaxation time and a(n) = I — 3e + Ae[(d x (f)) 4: + 
(9 z </>) 4 ] imposes a fourfold anisotropy with strength e in 
two dimensions. The function q((f) = (I — <fi)/2 dictates 
how the diffusivity varies across the interface. The inter- 
face thickness is given by wo while A is treated as numer- 
ical convergence parameter of the model. 

The parameters A, w a and t„ can be shown to be inter- 
related through the thin interface relations developed in 
Refs. |l2|, EH to map the above phase field model onto 
Eqs. CD)-© • Specifically, once a particular lambda is 
chosen, the thin interface relations establish a unique 
choice of Wq and To such as to yield the same d a and (3 in 
simulations. The aim is to choose a rather diffuse value 
of wo in order to expedite numerical efficiency In this 
work, we assume the interface kinetics coefficient ft w 0, 
to lowest order. This ability to quantitatively model the 
same materials parameters (3 and do is largely due to the 
term containing the constant a in Eq. ([5]). The term is 
called the so-called "antitrapping current", whose func- 
tion is to self-consistently counter the spurious effects 
of an interface thickness that is artificially enlarged for 
practical purposes. 



M (K/wt%) 


3.00 


CO (wt%) 


0.34 


k 


0.15 


D {^m 2 /s) 


3400 


T (K ■ ixm) 


0.10 


£ 


0.02 



TABLE I: Material parameters that define the samples em- 
ployed, m is the liquidus slope, Co the alloy composition, k 
the partition coefficient, D the diffusivity of impurities in the 
liquid, r the Gibbs-Thomson constant and e the anisotropy 
strength. 



The material parameters employed represent an Al-Cu 
alloy and are presented in tablefl] We solve the phase field 
equations using either a finite difference explicit Eulcr 
scheme on a uniform mesh, or, for simulations on larger 
scale, with a new finite difference adaptive mesh algo- 
rithm that utilizes a data structure developed Provatas 
and co-workers [IJ, [2l[ ■ 

Figure ([4]) shows an example of a typical sequence 
of directionally solidified dendrite arms growth under 
steady state conditions, i.e. a constant thermal gradi- 
ent (G — 5 K/mm) and puling speed (V p = 10 /zm/sec). 
In order to systematically study the evolution of interface 
structure, and, in particular, the selection of the colum- 
nar microstructure spacing, we also perform a power 
spectral analysis on simulated interface profiles using a 
Fast Fourier transform. An example of such evolution is 
also shown in figure (0}. It is noted that the emergence 
of a split in the main peak of the spectrum at the ear- 
lier stage reveals the interaction between dendrites that 
eventually leads to cell elimination. 

The phase field model was primarily used to simulate 
transient cooling conditions relevant to the experimental 
situation described in section (jTT|) . This was done by us- 
ing the thermocouple data to extract the local thermal 
gradient across the solid-liquid interface and the effective 
velocity of the solidification front, which were then fitted 
to provide functions used to determine G(t) and V p (t). 
The pulling speed was modeled after a fit of the front ve- 
locity obtained from experiments which, given that the 
interface is initially positioned at 7l, is systematically 
lower than the actual front velocity, the discrepancy de- 
creasing as the system evolves. Further details about this 
are discussed in section (IIVC|) . 

Figure ([5]) shows a typical spatiotemporal evolution 
of dendritic microstructures from our transient dynamic 
simulations. The morphological evolution has features in 
common with the experimental morphology in figure ^ . 
Specifically, at early times, when the velocity is fastest, 
the morphology is cellular and then starts to exhibit 
a series of kinetically-induced tip-splitting instabilities, 
giving rise to rising columns of seaweed-like structures. 
At later times, when the velocity decreases, the cellular 
branches emerge once more. The analysis of our phase 
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FIG. 4: (Top) Simulated evolution of a directionally solidified 
dendrite array grown under constant velocity and thermal 
gradient. Interface initially morphologically noisy. Distances 
are in pan. (Bottom) Spectral analysis of the last recorded 
interface structure (black line) as well as that of a close earlier 
time (red line). 



field simulations and the experimental data of section (|TT|) 
is discussed below. 



IV. RESULTS AND DISCUSSION 
A. Finite Size Effects 

In order to avoid boundary effects on the dendrite arm 
spacing, we study arm spacing in as large a single grain 
as possible. The simulations, on the other hand, are more 
time consuming for larger systems. To estimate a con- 
vinient grain (or system) size to use experimentally and 
theoretically in our spacing selection analysis, we studied 
the dependence of dendrite spacing on the system size in 
phase field simulations with constant control parameters. 
Figure ^ shows the final steady state spacing as a func- 
tion of transverse system size (i.e. grain size) for four 
different puling speeds. In all cases we started with a 
morphologically noisy interface. 

These results indicate boundary effects become negli- 
gible at systems (grain sizes) larger than about 10 3 /xra, 
even for the smallest rate of solidification, which leads 




FIG. 5: Dendritic microstructure evolved during the PF sim- 
ulation utilizing adaptive mesh refinement at the positions(a) 
600 (b) 1200 (c) 3600 and (d) 11400 microns from the initial 
position of the interface, corresponding to Tl- Each image 
is a zoom-in of the interface and is 65 microns in dimension. 
The colour represents concentration. 
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FIG. 6: Steady-state spacing versus transverse sample width 
for different pulling speeds. The effect of the boundary be- 
comes less prominent as the sample width increases. G=5 
K/mm. 



imental data does not compare well quantitatively and, 
especially, qualitatively with the unsteady state model of 
Hunt and Lu. Plausible reasons for this will be addressed 
in section (|V|) . 

Losert et al associate a rapid change or jump in spacing 
with the period doubling instability studied by Warren 
and Langer, citing boundary effects to account for the 
discrepancy of the jump being less than a factor of two. 
We propose that the discrepancy is also due to the fact 
that instabilities do not affect the system as a whole but 
are confined to smaller domains. The solute field ahead 
of the interface is not periodic and a localized distur- 
bance has to propagate along the front before affecting 
other regions. Hence there is a delay between an insta- 
bility in one region and its consequences to farther away 
dendrites. The spacing at any given time is thus the av- 
erage of spacing in different domains, rather than a con- 
sequence of an instability affecting the system instanta- 
neously. We examined this hypothesis in our experiments 
by examining different regions of a single transverse sec- 
tion at a specific position away from the chilling surface. 
Regions with different average spacing are found, consis- 
tent with our hypothesis. Figure <j8j) shows two separated 
regions within same transverse section of a grain which 
have different average spacing. 



to the largest spacing. The solidification rates studied 
experimentally and theoretically in the transient solidifi- 
cation analysis below are higher than those studied here, 
making finite size effects even more negligible in systems 
of 10 3 lira or larger. 



B. Dendrite spacing evolution: experiments 

Figure ([7]) shows experimental (black curve) and simu- 
lated (red curve) plots of the transverse dendrite spacing 
as a function of the front velocity, which varies during so- 
lidification. Also shown for comparison is the correspond- 
ing primary arm spacing determined from the unsteady 
state model of Hunt and Lu [23] (blue curve). In the 
experiments, temperature gradient is coupled with the 
velocity and are not independent. The solidification rate 
is obtained from the cooling curves by considering that 
a thermocouple registers a sharp change in the tempera- 
ture slope (with respect to time) when the solidification 
front passes through it. The experimental spacing shown 
is that obtained from transverse sections. Analysis of the 
longitudinal sections shows the same qualitative behav- 
ior, although the values are different, as expected. It is 
recalled that the experimental spacing data is obtained 
from different grains of the same experiment. 

The experimental results in figure are consistent 
with those reported by Losert et al |10| . Namely, the 
spacing exhibits regions of very slow change, between 
which it changes rapidly. It is noteworthy that the exper- 
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Velocity (mm/sec) 

FIG. 7: Spacing as a function of velocity obtained experi- 
mentally (black curve) and numerically (red curve). The blue 
curve corresponds to the relationship obtained by Hunt and 
Lu under unsteady state solidification. (Note that G(t) is also 
dynamically changing at each point.) 
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FIG. 8: separate zones within a grain at the same height level 
from the chilling wall which exhibit different values of average 
spacing. 



these tertiary branches as well as nucleation events. 




FIG. 9: Cast Al-Cu alloy showcasing the sharp change in 
morphology and spacing that occurs during a transition from 
low cooling rate (large-scale structure) to rapid cooling rate 
(finer structure). 



C. Dendrite spacing evolution: simulations 

In order to approximate the growth conditions simi- 
lar to those in our experiments, our directional solidifi- 
cation simulations included a variable thermal gradient 
and pulling speed, the values of which were set by fitting 
the corresponding curves obtained from experiments. Of- 
course, the pulling speed is not the same as the front 
velocity, especially under transient effects, since the in- 
terface moves within the coexistence region as growth 
conditions vary. This discrepancy is largest at the onset 
of the simulation due to the initial conditions. 

The two dimensional simulations in figure (O show 
remarkably similar behavior to the experiments, consis- 
tent with the idea that the branch spacings are confined 
to a band of slowly varying spacing separated by rapid 
change at particular velocities. The quantitative discrep- 
ancies between experiments and numerical results are to 
be expected due to the different growth conditions of the 
simulations compared to the experiments and that the 
simulations are two dimensional. 

It is noteworthy that the mean peak of the power spec- 
tra of the simulation data essentially captures the main 
branch spacing and not the seaweed-like sub-structure ev- 
idenced in figure J5]). The latter structures are present at 
early and intermediate times and likely arise due to inter- 
face kinetics induced by the high solidification rate. As 
mentioned above, these structures appear in the experi- 
mental data as well (see figure [2]). To illustrate this fur- 
ther, we cooled a sample at very low rate for a short pe- 
riod of time an then increased the cooling rate abruptly. 
The resulting microstructure is shown in figure ©. At 
low cooling rate, where the microstructure is larger, the 
microstructure comprises both cellular and seaweed-like 
tip splitting and side branching. After the cooling rate 
is increased, a finer structure emerges from the tips of 



V. MECHANISMS FOR SPACING EVOLUTION 

The initial pulling speed implies the diffusion length 
is several orders of magnitude smaller than the system 
width. Moreover, even though the pulling speed de- 
creases with time, the diffusion length stays several or- 
ders of magnitude below the system size. At the same 
time, the diffusion length is within an order of magni- 
tude from the dendritic spacing. Therefore, each den- 
drite branch can be considered to be in one of many, 
weakly interacting, sub-arrays of branches. The change 
in cooling conditions, on the other hand, is a global event 
that affects the whole system at the same time. It is 
therefore to be expected that the quasi-independent sub- 
arrays of branches will react differently to the globally 
changing cooling conditions. As a result, each sub-array 
may undergo instabilities -to some extent- independent 
from each other. This will influence the "jumps" in the 
spacing versus velocity curve. The rate of change in the 
cooling conditions may also play an important role, as 
the change in the global cooling conditions may cause 
the solute field associate with different regions to transfer 
information to neighbouring regions at different times. 

For a small number of dendrites branches in a system, 
the larger the change in the spacing due to a given branch 
elimination event. Conversely, a large enough system, 
or smaller the spacing, should exhibit a smoother evolu- 
tion of spacing. For the transient conditions we explored, 
there appears to exist quasi-independent domains which 
undergo instabilities without, for the most part, affecting 
each other. After an instability, such as cell elimination, 
affects one, or several, entire domains, the average inter- 
dendritic spacing of the entire system will experience a 
sharp change. On the other hand, we also expect pe- 
riods of low latency in which the overall spacing does 



9 



not change, or evolves very slowly, while the different re- 
gions slowly interact through diffusion of the solute field. 
These observations are consistent with our findings and 
also with those of Losert et al [lOj , in which he changed 
the pulling speed in steps and found a sharp jump in the 
spacing after it went almost unchanged for a range of 
pulling speeds. 

The existence of plateaus connected by rapid changes 
in spacing points to a dynamics in which an energy bar- 
rier has to be overcome for the spacing to adapt. Specif- 
ically, for the tip of a cell or branch to split, it effectively 
needs to pass through a flattening stage in which the 
tip radius becomes effectively infinite. This lowers the 
interface undercooling (i.e. the contribution of dendrite 
tip curvature to interface undercooling, doK , becomes 
zero), preventing tip splitting until a larger local inter- 
face velocity is reached. Increasing the velocity reduces 
the diffusion length and shortens the distance over which 
solute is rejected. This effect acts to reduce arm spacing. 
Hence, there will be a competition between the two ef- 
fects, keeping the spacing constant due to the first effect 
until the driving force provided by second phenomenon 
is large enough to force dendrites or cells to split and re- 
duce their spacing. Figure (|10[) schematically represents 
the stages of tip splitting. 
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FIG. 10: Schematic representation of different stages during 
the tip splitting mechanism. In the flattening stage (b), the 
undercooling due to the Gibbs Thomson effect is essentially 
zero. 

An increase in branch spacing via cell elimination oc- 
curs when a dendrite is blocked by either the secondary 
arms (in low speed solidification) of an adjacent dendrite 
or by the main stalk (in high speed solidification) of an 
adjacent dendrite. As shown in figure (fTTj) . in both cases, 
the seizing arm grows at an angle with a velocity compo- 
nent perpendicular to the growth direction of the primary 
arms (which is the direction of lowest energy, obtained 
as the product of anisotropy direction and heat flow di- 
rection). Unlike equiaxed growth, where secondary arms 
grow at a rate comparable to the growth rate of pri- 
mary arms, in directional solidification the growth rate 
of secondary branches is negligible compared to that of 
the main trunks. Thus, the seizing mechanism described 



above is not able to act until diffusion-mediated inter- 
actions through the melt become significant. As growth 
velocity decreases, solute diffuses a longer distance and 
if the diffusion length of a dendrite is large enough, it 
interacts with adjacent diffusion fields. It is at this point 
that a jump in system energy can occur. 




fa) 

FIG. 11: Cell elimination caused by the neighbouring (a) pri- 
mary arms or (b) secondary arms [231 ] at high and low velocity 
solidification respectively. 



Hence, in the case of cell elimination, there is an en- 
ergy increase due to the growth in any of the undesired 
directions and an energy decrease for changes that in- 
crease the distance between dendrites (i.e. characterized 
by non-interacting diffusion fields). These two phenom- 
ena dynamically compete, with the spacing change even- 
tually determined by the dominant effect, where the for- 
mer effect acts as a barrier against the change in the 
spacing and the latter provides the driving force for cell 
elimination. It is also plausible that over a small range 
of cooling conditions (i.e. solidification rate, temperature 
gradient) these effects may balance each other, causing 
the spacing not to change very much, at least over some 
long-lived transient time (i.e. the plateau regions in the 
data). Furthermore, the farther from steady state spac- 
ing the initial condition of the system, the higher the 
energy and the larger the driving force required for the 
system to perform branch elimination. Thus, it is pos- 
sible that different spacings can also exist at the same 
cooling conditions depending on the history of the sys- 
tem. This matter will be further investigated in future 
work. 
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